Effect of Subband Landau Level Coupling to the Linearly Dispersing Collective Mode 

in a Quantum Hall Ferromagnet 
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In a recent experiment (Phys. Rev. Lett. 87, 036803 (2001)), Spielman et al observed a linearly 
dispersing collective mode in quantum Hall ferromagnet. While it qualitatively agrees with the 
Goldstone mode dispersion at small wave vector, the experimental mode velocity is slower than 
that calculated by previous theories by a factor of about 0.55. A quantitative correction may be 
achieved by taking the subband Landau level coupling into account due to the finiteness of the 
layer thickness, which gives a better agreement with the experimental data. A method coupling the 
quantum fluctuation to the tunneling is briefly discussed. 



PACS numbers: 73.43.-1,71. 35.Lk, 73. 21.-b, 73.40. Gk 

The fractional quantum Hall effects as well as most 
of related phenomena reflect the behavior of the ground 
states of the two-dimensional electron gas in high mag- 
netic field at zero temperature ||l|,||. However, a first ex- 
perimental signal of the finite temperature phase transi- 
tion in such systems has been recently observed by Spiel- 
man et al in a bilayer system for i^t = 1 [HI]- This 
transition seems to be closely related to a Kosterlitz- 
Thouless transition into an earlier predicted broken sym- 
metry state Q which was akin to the Josephson tun- 
neling in superconductivity. However, differing from the 
Josephson effect in which the zero-bias conductance at 
the zero temperature was divergent, the conductance 
peak in that experiment was finite even extrapolating 
to the zero-temperature. This phenomenon, with other 
intriguing novel properties, caused a set of theoretical re- 
search works §,g,|,|jl^,|l|,|l|,|l|,|l|,^ in which various 
scenarios have been suggested while there is still a variety 
of open issues ( for a short review, see p6|). 

Among these issues, we shall focus on the finite layer 
thickness affecting the tunneling through the tilted field, 
which was disregarded by literature. As a result from the 
zero thickness approximation, a quantitative discrepancy 
between the previous theoretical calculation and the ex- 
perimental data already appeared. The precise descrip- 
tion of the issue is as follows: In a recent experiment [Q, 
Spielman et al observed a linearly dispersing collective 
mode in a small wave vector in bilayer two-dimensional 
electron systems, which was identified as the pseudospin 
Goldstone mode long expected. However, Figure 4 of 
their published paper showed the sound velocity, namely, 
the slope of the linear dispersion experimentally observed 
was smaller than that theoretically predicted by a factor 
of about 0.55. One argued that this discrepancy may re- 
sult from the overestimate in the Hartree-Fock approxi- 
mation used in the theoretical calculation, and the quan- 
tum fluctuation may repair this discrepancy. However, 
the present exact diagonalization result indicated that 
the correction from the quantum fluctuation does not 
seems to reach such a small factor Thereby, it is 

worth to look for other sources to influence the mode. 



On the other hand, due to the dispersion proportional 
to the in-plane magnetic field, the finite thickness of 
the layers may affect it because the layer thickness is 
comparable to the interlayer spacing in the experiment. 
However, the previous theoretical treatments were essen- 
tially based on the zero-layer- thickness approximation. 
The experiences in the study of the quantum Hall sys- 
tems were that finite thickness effects often determine the 
quantitative consistence between the theory and experi- 
ment, e.g., the gap of Laughlin's state, and so on. 

Moreover, for this bilayer system, the ground-state be- 
havior is not well understood yet if there is a tilted field 
and the finite layer thickness is taken into account. We 
even do not have a satisfactory variational ground-state 
wave function. Thus, one may suspect whether this ex- 
perimental observed mode can be identified as the theo- 
retically anticipated pseudospin Goldstone mode. In this 
paper, we deal with the experimental data through a the- 
oretical model which attributes the finite thickness cor- 
rection to the sound velocity. If the correction is positive, 
we may say that the identification can be accepted and 
otherwise it could be more suspect. In a composite boson 
picture when the field is tilted |jl^, while the composite 
bosons see an opposite effective perpendicular field with 
an equal magnitude in the mean- field state, they see 
a weakened effective parallel field B*^ due to the finite 
layer thickness via the subband Landau-level coupling. 
Thus, one can apply the the interlayer tunneling theory 
P,p]pX|] to this composite boson system but the parallel 
field _B|| is reduced to B*^ By using the Hartree-Fock 
estimate of the sound velocity jl^JT^, one can obtain 
the linear dispersing Goldstone mode in the same way 
as the theoretical line in Fig. 4 of Q but q is replaced 
by q* = eB^^d/h. This dispersion has a substantial im- 
provement to fit the experimental data. Thus, we have 
a positive correction. In the meanwhile, our composite 
boson formalism shows there is a coupling of the quan- 
tum fluctuation to the tunneling. However, a detailed 
discussion of the quantum fluctuation has exceeded the 
goal of the present work. 

In this work, we only deal with the finite-thickness ef- 
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feet via the subband Landau-level coupling and neglect 
the others, say, affecting the Coulomb interaction due to 
the finite thickness. In order to deal with the subband 
Landau-level coupling analytically, we assume that the 
electron gas is confined in a plane by an infinite harmonic 
potential. Before going to the details, we argue that this 
choice of the confining potential may quantitatively cor- 
rect the theoretical calculation without severely impact- 
ing the comparability between our calculation and the 
experimental data despite the fact that the realistic con- 
fining potential in the sample is essentially a finite square 
well. The harmonic well is very different from the square 
well in their excited spectra, for the harmonic spectrum 
is equal gaped while that of the square well is not. How- 
ever, the temperature is extremely low so that only the 
lowest subband of the lowest Landau level is filled in the 
present situation. Thus, no excited spectra will be in- 
volved. In the ground state, one may variationally adjust 
the harmonic frequency such that the ground state wave 
function has the best shape, and let the subband energy 
equal that in a more realistic square well. In this sense, 
the harmonic potential may be a good approximation to 
a realistic potential, to give a quantitative correction due 
to the subband coupling. Such a harmonic potential has 
been chosen to deal with many quantum Hall systems to 
replace the realistic potential which is either triangular 
pof or square It was also used to discuss a giant 

magnetoresistance induced by a parallel magnetic field 

We start from the problem of a single particle in 
a strong magnetic field which is tilted at an angle to 
the x-y plane. An in-plane field in the x-direction vio- 
lates the two-dimensional rotational symmetry. By in- 
troducing a harmonic confining potential with the char- 
acter frequency O in the z-direction, the electron is 
restricted to quasi-two-dimensions. The single-particle 
Hamiltonian can be diagonalized as Hg.p. — huj^a^^a^ + 
fiuj+alaz with the diagonalized oscillators given by — 
(a^, al, a^, az)X~'^ , a\ = (a^, o|, a^, az)X^'^ , where = 
-^{-d, + m^a, = ^(%+ K) and 4 = ^(-J. + 
z'),az = + z'). z' = n^/'^z. The vectors X± are 

given by 



X CX {lOc + UJ-, - 



X+ oc (- 



UJc + 1-1-'+ 



The frequencies uj± are given by uj\ — ^(17^ 4- w^) ± 
\\J {9? - ^2)2 + 4|(I;pl7wc, where Cb = Wa;(wc/i7)i/2 g^^^ 
172 = 172 + tj2. g^j^g -j-j^g cyclotron frequencies 

corresponding to Bx and Bz- Here we have applied the 
unit Ib — y/Tic/eBz = 1. In addition, there is a conser- 



vation quantity — aj^a^ with cll — "^(^C + ^-iid 
\ = -^{—d^ + irrespective of whether the tilted 



angle 9 = or not. To solve this single-particle prob- 
lem, we seek the ground state which is the eigenfunction 
of L^. It is useful to make a coordinate rotation with 
^^^ = £,+P^+ -fz' and z' — z' with (3 and 7 deter- 
mined by [a^,£,] = [a^,^] = 0: 
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4 J^3 



{X^ - X^)X^ - {X2 - X 
{X2 - X4)X-^ - {X2 - )X^ ' 
2(^2 X-^ — X^ X-^ ) 

(x+ - xt)x- - (x- - x-)x+ ■ 



(2) 



The ground state wave functions are highly degenerate 
and of the form ^o{i,t,z') = /(|)e9 with g{lC,S') 
being a quadratic form of ^,$*, z' whose coefficients are 
determined by a^e^ = a^e^ = 0. The function /(^) is an 
arbitrary function of ^. 

Notice that linear-independent wave functions ^™e^ 
(m=0,l,2,...) are not the eigen functions of L^. How- 
ever, one can start from those linear-independent wave 
functions to construct the common eigen functions of 
Hs.p. and Lj, which read fm{C)e^i with /m(0 = 
X)m'=o fmm'S,"^ for M being the number of Landau orbits 
and M = N for j/y = 1 . The coefficients fmm' are depen- 
dent on the in-plane field and confined by /mm(0) = 1 
and /mm'{0) = for TO ^ to' if = 0. Those degenerate 
ground state wave functions are orthogonal and with the 
eigen value to of L^. 

After solving the single-particle problem, we turn to 
the many-body ground state wave function. The com- 
mon Laughlin's state or Halperin's (lll)-state for the 
vanishing tilted angle is no longer a good variational wave 
function. However, to be enlightened by them, we pos- 
tulate the many-body ground states for vt — ^ SiS 

^-0(7^1,...,^^) = St(/o,.../jv-i)exp(^.g,) x |PSS), 



where are the three-dimensional position vectors and 
St is the Slater determinant of /o('Cp(i)), fN-iiS,p{N))- 
|PSS) denotes the possible layer coherence, which is ei- 
ther pseudospin fully polarized or coherent. Although the 
quantity L = L^. is no longer conserved due to the in- 
teraction, there is a conserved quantity L + • 
The constructed state is the eigenstate of this quantity 
with eigen value 0. Due to Af = N, one has exactly 

vE'o(fi,...,fAr) CX Hte -0)exp(^5.) X |PSS), (3) 

i<j i 

where, for the balanced bilayer system, i — l,...,N/2 
and 1 -|- N/2, ...,N denote the particles in layers 1 and 2, 
respectively. The spatial part of this wave function goes 
back to Halperin's (lll)-state for = 0. Corresponding 
to the system in the experiment of Spielman et al 
the coherent state is given by 



|PSS) = 



(4) 



>) = ^(|T)+e^^|i)), 
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if one uses | tii) to represent the electron in the upper 
or lower layer. 

Since the electrons in the quantum Hall state are 
strongly correlated, when an electron tunnels from one 
layer to another, the fluxes combined with the electron 
also move accompanying the electron. If the magnetic 
field is not tilted, the fluxes accompanying the electron 
do not move in the tunneling because they are perpendic- 
ular to the x-y plane. However, if the field is tilted, the 
fluxes accompanying the electron no longer lie in the x-y 
plane due to ^ ^ ^. Thus, the electron tunneling causes 
the flux hopping from one layer to another. This can not 
be reflected in the single electron tunneling picture. A 
better formalism is the composite boson formalism. The 
single composite boson tunneling counts the charge tun- 
neling and the flux hopping simultaneously. Hence, we 
use the composite boson formalism. According to @ , for 
a, vt — ^ state, the composite boson theory in a tilted 
magnetic field can be achieved by the anyon transforma- 
tion in the spatial wave function [[L8[ 



Kij I 



(5) 



In terms of Ref. [Q, this transformation gives a sta- 
tistical gauge field, af_,{n) = -J2j^t-ff-p-^0'z{n) = 

Ej^j lIT^' where i = S: + iy; ii2 = I + f3 and 
£21 = — 1 + /3. The corresponding statistical magnetic 
field 6 = V X a, i.e., 

=-27r(l + /3)-i(f]K)V2^^<5(2)te,). (6) 

In the mean field approximation, — 2Bz, while 6|| 
partially cancels -By . The effective parallel magnetic field 
seen by the composite boson reads 



(7) 



Bl = {l- (7/1 + l3)\J^ujc/uj\\)B\\ = DB 



We shall see that D <1 and the equality holds liuOc/Vl 

3 

0. Since 7/(1-)- /?) cx ^ + 0(^) for < Wc, £> is almost 
independent of _B|| for the parameters used in The 
inter-layer tunneling operator reads 

T = -t J dxdy^p^^{x,y)e''^''^^pi{x,y) + h.c. 

= -t I dxdycf>\{x,y)e'^^'--'')Ul{x,y) + h.c., (8) 



where ip^ and denote the electron and composite boson 
creation operators in the upper layer respectively, etc. d 
is the interlayer spacing. Az is the the 2:-axial component 
of the vector potential of the external magnetic field and 



Qz is that corresponding to b{r). Thus, in the mean field 
approximation, 4>t,i = y/poe^^^-^ = —^i), Az — Uz = 
DBi\x and the tunneling Hamiltonian reads 



T = — I dxdy— cos((y5 



q X) 



(9) 



where the effective wave vector q* = Dq instead of 
q — 27ri?||(i/$o- Thus, the theory |0 applied to the elec- 
tron can also apply to the composite boson because the 
anyon transformation is carried out in the spatial wave 
function and does not change the particle density and the 
pscudospin coherence state. The only difference in the fi- 
nal result from that in the electron's case is replacing q 
by q* . That is, for small q* , the linear dispersion is given 
by vq* . However, all measurements in the experiment 
by Spielman et al Q correspond to namely, q. To 
compare with the experimental data, the dispersion may 
written as 



vq* — Dvq ~ v*q. 



(10) 



Due to D < I, V* is smaller than u by a factor D. If 
we use the Hartree-Fock value of the sound velocity to 
estimate v, the dispersion is exactly what is plotted in 
Fig. 4 of Q except that q is replaced by g*; or in the 
wave vector q, the sound velocity v is renormalized non- 
perturbatively to v* . This renormalization is remarkably 
different from that by the quantum fluctuation in the 
zero-thickness theory and is the result of the strong cor- 
relation of the electrons. 

It is easy to see that the magnitude of D is determined 
by the ratios wy/wc and lOc/^- Since the tunneling peak 
is destroyed in a small tilted angle, say ~ 0.6 T (thus 
the critical value of the ratio is wy/wc ~ 0.3 for ut = 
5.2 X lO^^/cm^ because Bz 2T ) one can consider 
the small ratio Wy/Wc only. In fact, one finds that D is 
almost independent of B\\ if lu\\/uJc < 0.2. In Table 1, we 
list the magnitudes of D for uj\\/uJc = 0.2. 



UJc 

0. 


0.10 


0.40 


0.50 


0.80 


1.00 


1.25 


1.67 2.50 


D 


1.00 


0.90 


0.86 


0.76 


0.70 


0.65 


0.60 0.55 



Table 1 The magnitudes of D for the different ratios ^ 
and ^ = 0.2. 

The first column in Table 1 implies that it is back 

to the result got by ignoring the finite thickness. The 

last value of D reaches the experimental data but the 

corresponding ratio jf = 2.5 is too large for the sample. 

In the experiment sample, the square well has the width 

a = ISnm. The magnetic length Ib ~ 17nm for the 

density ut — 5.2 x 10^° and i^t — 1- The ratio fl/ojc 

is dependent on the heights of the square well. For the 

potentials with the infinite height on two sides of the well, 
2 

ujc/^ ^ -§75- ^ 0.11 and D — 0.99 if we determine V, by 
If this was the case, there would 



En = 



2 ^' 2ma^ 

be almost no finite tliickness effect. However, the wells 
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in the real bilayer systems are finite and asymmetric. 
The single particle ground state energy may be lowered 
substantially. If the potential Vi on one side is lower than 
V2 on the other side, the ratio is given by 

where k < ki is determined by the solution of the equa- 
tion ka = n — sin~^(A:/A:i) — sin~"'^(fc/fc2) with ki — 
y^2mVi/h. Thus, the ratio is dependent on the heights 
of the well on both sides. Assuming Vi = V2 and for a 
typical well with Vi ~ /{ma?), the value of k^ may 
reduce to 0.157r^/a^. Hence, 

uJc/VL^a^/{0.l^'KH]j)^Q.7 and D - 0.78. 

The details of the well heights were not reported in . 
We do not expect that the result stemming from the 
mean field approximation can completely fit with the ex- 
periment. However, what we can say is that the finite 
thickness correction drives the theoretical results to a 
positive direction to be consistent with the experimental 
data. There may be other sources to improve the theo- 
retical result. Especially, the quantum fluctuations may 
be very important because it lowers the spin stiffness and 
then flatten the slope. Typically, the quantum fluctua- 
tion may lower the stiffness by a factor < 0.80 and then a 
factor < 0.89 to the sound velocity according to an exact 
diagonalization result by Moon et al |l7| . Combining this 
factor with the factor from the finite thickness correction, 
the theoretical result may have a good agreement with 
the experimental data. 

Before finishing this paper, we would like to point out 
that the composite boson formalism we used here intro- 
duces a coupling of the phase to the quantum fluctu- 
ation other than that discussed in literature jl7|. Note 
that from (||) to (^), we use the mean field approximation. 
That is, 5a z in ^ — Az — 5a z has been neglected. 
The gauge fluctuation couples to the phase in the tunnel 
Hamiltonian is given by 



r = -— [cos(^ 

ZTT 



q* x) cos{5azd) 
sin((p — q*x) siii{5azd)]. 



(11) 



Obviously, the gauge fluctuation will weaken the tunnel- 
ing, then the charged gap AgAS- We will leave further 
discussion of this aspect to a separate work . 

In conclusion, we have considered the subband Landau 
level coupling in the interlayer tunneling of the bilayer 
quantum Hall system. It is found that this coupling may 
quantitatively affect the linear dispersion of the pseu- 
dospin Goldstone mode. A further improvement to the 
linear dispersion requires more sample parameters, espe- 
cially, the heights of the confining potential as well as to 
develop a reliable method to count the correction from 
the quantum fluctuation. 

This work was supported in part by the NSF of China. 
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